The Flux-Flux Correlation Function for Anharmonic Barriers 
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The flux-flux correlation function formalism is a standard and widely used approach for the 
computation of reaction rates. In this paper we introduce a method to compute the classical and 
quantum flux-flux correlation functions for anharmonic barriers essentially analytically through the 
use of the classical and quantum normal forms. In the quantum case we show that the quantum 
normal form reduces the computation of the flux-flux correlation function to that of an effective 
one dimensional anharmonic barrier. The example of the computation of the quantum flux-flux 
correlation function for a fourth order anharmonic barrier is worked out in detail, and we present 
an analytical expression for the quantum mechanical microcanonical flux-flux correlation function. 
We then give a discussion of the short-time and harmonic limits. 
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I. INTRODUCTION 

In spite of the tremendous increase in computer power 
over the last decades, the computation of quantum re- 
action rates still is a formidable task. This is the topic 
of this paper, and we begin by giving the setting that is 
relevant to our work. 

The microcanonical rate constant is given by 



k{E) = (2Trhp r {E)) 1 N(E), 



(1) 



where p T (E) is the density of states of reactants and 
N(E) is the cumulative reaction probability, which in 
turn can be formally expressed in terms of the S matrix 



(2) 



Here the inner sum runs over the asymptotic states of 
reactants and products which are labeled by n T and n p , 
respectively, and the outer sum covers all values of total 
angular momentum J. 

Though formally correct, the computation of a reac- 
tion rate via the S matrix is extremely inefficient since 
the computationally expensive information related to the 
state-to-state reactivities embodied in the S matrix is 
" thrown away" as a consequence of the averaging embod- 
ied in the summations in $Z§. Motivated by the success 
of transition state theory (TST) for computing reaction 
rates using classical mechanics, many researchers have 
sought a quantum mechanical version of transition state 
theory. Recall that the main idea of TST as invented by 
Eyring, Polanyi and Wigner in the 1930s is to compute 
classical reaction rates from the flux through a dividing 
surface in phase space which separates the phase space 
region associated with reactants from the phase space 
region associated with products. Assuming the dividing 
surface to be given by an equation s(z) = 0, where s 
is a scalar function on the 2 /-dimensional phase space 



with coordinates z = (q, p) = (qi, . . . , qf,pi, . ■ ■ ,Pf) 
with s(z) < in the reactants region and s(z) > in the 
products region, the classical microcanonical rate con- 
stant can be written as 

k CL (E)^ Pl (E)- 1 (2Trh)-f I dz S(E - H(z))F(z) . (3) 



Here F(z) is the so called flux factor which is given by 



F(z) 



dt 



= -{H,Q(s)}(z) (4) 



with denoting the Heaviside function, {•,■} denot- 
ing the Poisson bracket and z t denoting the solution of 
Hamilton's equations at time t with initial conditions 
zq = z, i.e. Zj = $^(z), where <&h is the Hamiltonian 
flow which acts on z for time t. 

For the TST computation of the rate constant to be 
useful the dividing surface needs to have the property 
that it is crossed exactly once by reactive trajectories 
(i.e. trajectories evolving from reactants to products) 
and not crossed at all by nonreactive trajectories. A di- 
viding surface not satisfying this no-recrossing property 
leads to an overcstimation of the reaction rate. The con- 
struction of a recrossing free dividing surface has posed 
a major problem in the development of TST. Formally, 
recrossing trajectories are eliminated by multiplying the 
integrand in ^ by the projection function 



P r (z)= lim 9(a(z t )) 

t— >-\-oc x ' 



(5) 



which evaluates to 1 if the trajectory z t = $^ (z) evolves 
to products for t —> +oo, and to otherwise. It is then 
not difficult to see that the this way corrected expression 
for the rate constant can be written in the form 



/oo 
dtCcL (£,*), 
-oo 



(6) 
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where Cch(E,t) is the flux-flux correlation function 
(FFCF) 

Ccl{E, t) = (2nh)-f J dz 5(E - H{z))F(z)F(z t ) . (7) 

The canonical analogues of the microcanonical expres- 
sion above are easily obtained from replacing the den- 
sity of states S(E — H(z)) by its canonical counterpart 
exp(— pH(z)) (j3 = (k-QT)^ 1 denoting the inverse tem- 
perature). Miller, Schwartz and Tromp [l[ and Ya- 
mamoto Q took this as a starting point to express the 
quantum mechanical rate constant in a similar fashion. 
For the microcanonical case, this amounts to replacing 
the phase space functions above by the corresponding op- 
erators and the integrals over phase space by the traces 
of operators. For the quantum analogue of the rate con- 
stant in ([3]), this leads to 



k QM (E) = p T (E)~ lr Tr {d(E - H)FP r } , 



(8) 



where p r (E) now is the quantum mechanical partition 
function of the reactants (for which we use the same sym- 
bol as in the classical case), the flux factor (Q| becomes 
the operator 



dt \ 



t=o 



= -[H,e(s)] (9) 



and the projection function ([5]) becomes the operator 



Pr 



lim e 

t— >+oo 



s c 



-iHt/ft 



(10) 



Similarly to the classical case expression (jHJ for the quan- 
tum rate constant can then also be rewritten in terms of 
a flux-flux correlation function, namely 



k QM (E) = (2 Pr (E)Y 



dtC QM (E,t) (11) 



with 



Cqm(E, *) = Tr U( E - H)Fe i6t/h Fe- i6t / h \ . (12) 

Clearly, the ability to compute the quantum mechani- 
cal flux-flux correlation function is essential for comput- 
ing quantum mechanical reaction rates using this ap- 
proach. Exact expressions for the canonical quantum 
mechanical flux-flux correlation function can be found in 
[l[ for the free particle and in [l], [§, ,4] for the one degree- 
of-freedom (quadratic) parabolic barrier. Calculations 
for other one degree-of-frecdom systems have also been 
carried out numerically. In particular, the quantum me- 
chanical flux-flux correlation functions for the symmet- 
ric and asymmetric Eckart and the double- well potential 
have been studied in detail (see d, @)- However, these 
computations were no ab initio computations leading to 
analytical expressions, but were carried out using various 
computational techniques. Indeed, the development of 



computational techniques to compute the quantum me- 
chanical flux-flux correlation function has been a subject 
of great interest since the development of this approach. 

The traditional approach of quantum mechanics in- 
volves basis function techniques, and a review of this 
approach is given in 0, H[ . While this particular compu- 
tational approach has proven successful for "moderately" 
sized molecules, difficulties are encountered when "large" 
systems (such as biomolecules) are considered. The ob- 
stacle has been termed the "exponential wall" of difficulty 
that one encounters when attempting to perform numer- 
ical quantum calculations in the traditional manner for 
many degree of freedom systems 0] . Scmiclassical meth- 
ods provide an alternative approach. They have been pi- 
oneered by Miller and provide significant computational 
advantages for "large" systems; see [l(| for a recent re- 
view. 

In this paper we present a third approach based on 
recent work on the phase space approach to quantum 
mechanics that relies heavily on the dynamical systems 
framework arising from the application of classical and 
quantum normal form theory. Classically, this approach 
has provided a way of constructing a recrossing free divid- 
ing surface in the neighborhood of a saddle type equilib- 
rium point which forms the barrier between reactants and 
products (more precisely the existence of the recrossing 
free dividing surface is guaranteed for energies not too far 
away from the energy of the saddle) [HQ. 

The dividing 

surface at a given energy is bounded by a so called nor- 
mally hyperbolic invariant manifold (NHIM) [l5j which 
as an invariant subsystem with the given energy forms 
a 'supermoleculc' localized between reactants and prod- 
ucts which can be viewed as the transition state that gave 
TST its name. The dividing surface and the NHIM can 
be constructed in an algorithmic fashion based on a clas- 
sical normal form (CNF) which can be computed in the 
neighborhood of the saddle. This is a nonlinear canonical 
transformation of the phase space coordinates which, to 
any desired order, decouples the dynamics near the sad- 
dle into a single reactive mode and bath modes. (Note: 
the "decoupling" occurs through a local integrable ap- 
proximation valid in a neighborhood of the saddle, and 
does in general not imply separability.) More recently 
it has been shown that the classical construction can be 
generalized to the quantum case [l6|-[l8j . This essentially 
involves a systematic quantization of the canonical trans- 
formation leading to the decoupling in the classical case. 
The result is a so called quantum normal form (QNF) 
which leads to a local decoupling also of the quantum 
dynamics into a reactive mode and bath modes. 

The main objective of this paper is to use the classi- 
cal and quantum normal forms to analytically compute 
the classical and quantum FFCFs for anharmonic barri- 
ers. The outline is as follows. In Sec.|H]we describe how 
the classical and quantum normal forms can be used to 
express the classical and quantum mechanical FFCFs. 
Sec. HI] A is concerned with the classical case, and we 
show how the dividing surface having the no-recrossing 
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property obtained from classical normal form theory en- 
ables a computation of the FFCF that does not require 
the computation of trajectories. Sec. [II] B is concerned 
with the quantum case and we show how the quantum 
normal form approach enables a reduction of the com- 
putation of the FFCF in the general / degree-of-freedom 
case to the case of a one degree-of-freedom anharmonic 
barrier. Then in Sec. lIIII we work out in detail the micro- 
canonical quantum FFCF for the example of a barrier in 
one dimension described by a fourth order anharmonic 
barrier. 

By itself this result should be of some interest because 
it is the first explicit example of an exact calculation 
of the FFCF for an anharmonic potential. In Sec. IIVI 
we give a discussion of our results and an outlook for 
future work. The technicalities of two computations are 
deferred to appendices. 



II. THE FLUX-FLUX CORRELATION 
FUNCTION IN THE FRAMEWORK OF 
CLASSICAL AND QUANTUM NORMAL FORM 
THEORY 

We begin by showing how classical and quantum nor- 
mal form theory can be used to calculate both the clas- 
sical and quantum FFCF for reactions associated with 
barrier given by an index one saddle of the potential en- 
ergy surface. In fact the (slightly more general) start- 
ing point is an equilibrium point of Hamilton's equations 
which is of saddle-center-.... -center stability type. For a 
system with / degrees-of-freedom, this means that the 
matrix associated with the linearization of Hamilton's 
equations about the equilibrium point has one pair of 
real eigenvalues ±A and / — 1 complex conjugate pairs 
of imaginary eigenvalues ±iu>fc, k = 2,...,/. We will 
call such an equilibrium point a saddle for short. For 
simplicity, we will restrict ourselves to the generic situ- 
ation where the linear frequencies uik are not resonant, 
i.e. 7712W1 + . . . + mtuif ^ for every nonzero vector of 
integers (7712, . ..,m/). 



A. The classical case 

Classical normal form (CNF) theory [H Q3, QJ] pro- 
vides an algorithm for constructing a (nonlinear) canon- 
ical transformation z M> Z = (Qi, . . . , Qf, Pi, . . . , Pf) 
which, after truncation at a suitable order TV, leads to 
an integrable approximation of the dynamics near the 
saddle. In terms of the normal form coordinates Z the 
iVth order CNF of the original Hamilton function H(z) 
assumes the following form 

#cnf( z ) = Kcnf(I> hi Ja,---, Jf) 

LA72J fl ox 
n=0 \a\=n 



Here ['J denotes the floor function, the a = (ai, . . . , a/) 
are vectors with nonnegative integer components and 
norm \a\ = Y/ k =i a k, 



{ p i - Qi) 



(14) 



is an action type integral associated with the reactive 
mode, and 



Jk 



Ql) 



2,...,/, 



(15) 



are action integrals of the bath modes which we group 
together in the vector J. The CNF transformation (in- 
cluding the coefficients K n a in (|13[) ) can be computed in 
an algorithmic fashion as described in detail in (l7| . 

In terms of the normal form coordinates the dividing 
surface is given by s(Z) = Q\. Then, following the defi- 
nition in Eq. ((4]) the flux factor is given by 



F(Z) = -9(0, 



5{Q 1 )u(I,3)P 1 , 



where 



(16) 



(17) 



Following the FFCF then takes the form 
C C L(E,t) = [ dZ6(E-H { c ^ F (Z)) 



x5(Q 1 )S(Q lt )^ 2 (I,3)P 1 P lt . (18) 

The product of the (5-functions of Qi and the corre- 
sponding time evolved coordinate Q\ t (using the flow 

associated with -Hcnf) indicates that only the infinites- 
imally short time scales t — > give a non-vanishing 
contribution to the integral. The short-time expansion 
Qit = Qi + v(I,J)Pit + 0{t 2 ) and P u = Pi + 0(t) 
yields 



CcL(E,t) = 2f(E)8(t) 



(19) 



where 



= (27T)''- 1 f dIdJS(E - K { C N NF (I,J)MI,J) 

JR f . 



(27T)'- 1 / dJ, 
Jl(E,J)>0 



(20) 



with I = I(E, J) solving the energy equation 



#£nf( J > j ) = E. The last integral in Eq. {2D) is the vol- 
ume in the space of the center actions J = (J2, . . . , J/) 

enclosed by the contour /f^p(0,J) = E, and accord- 
ingly f(E) is nothing but the directional flux through 
the dividing surface [jji HH • 
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We note that formally the result given by Eq. (|T9|) ex- 
ists in the literature before (e.g., see Refs. d,[l^{ for a cor- 
responding canonical version of the formula). However, 
the contribution of the classical normal form theory in 
providing a dividing surface with the no-recrossing prop- 
erty provides a new formula, and interpretation, of the 
pre- factor f(E) in terms of an integral, in the bath mode 
action space, over the NHIM. This eliminates the need 
to compute trajectories (and the projection function (|SJ)) 
in the computation of the classical FFCF. 



B. The quantum mechanical case 



In [16l - ll8| a quantum normal form (QNF) procedure 
has been developed that yields a local decoupling of a 
reactive mode and the bath modes also in the quantum 
mechanical case if the corresponding classical system has 
a saddle equilibrium of the form described above. In 
the quantum case the local simplification of the Hamil- 
ton operator is achieved by conjugating it with a suit- 
able unitary transformation. Similar to the classical case 
this unitary transformation and the transformed Hamil- 
ton operator can be computed in an algorithmic fashion. 
The transformed operator then takes the form of a power 
series in terms of elementary operators associated with 
the reactive and bath modes and, in addition, Planck's 
constant. Truncating this expansion at a suitable order 
N gives the Nth order QNF approximation Hqn F which 
is of the form 



H, 



(N) 



~ ^QNf(^' ^2, J$, ■ ■ ■ , Jf) 

L7V/2J 



— ^ ^ l^n.ot ,j I ^ J*2 ' ■ ' ^ f ' (21) 

rc=0 \ a \+j=n 

Here, the notation is the same as in (|13p. where in addi- 
tion the j are nonnegative integers and 



i = \[ 



(22) 



is an operator associated with the reactive mode, and 



Jk = l[P£ + Q 



k = 2,...,f 



(23) 



are operators associated with the bath modes. In (f2"2")l 
and (|2U)) the Qk and Pk, k = are as usual pairs 

of conjugate position and momentum operators that sat- 
isfy the commutation relations [Qk,Qi] = [Pk->Pl] = 
and [Qk,Pi\ = ihSki- The approximation of the orig- 
inal Hamilton operator by the QNF in Eq. (f2~Tj) holds 
locally in the vicinity of the saddle equilibrium of the 
correspon ding classical system in a sense that is made 
precise in [171 ]. 

Since the trace of an operator is invariant under 
unitary conjugations of the operator we can evaluate 



Eq. (Ull) using the QNF to get 

Cqm(E, t) 

^'QNF I ^QNF 



Tr {S(E- e*SW* e^^} 

(24) 



with the flux operator given by 



F QNF = ^[^QNF> (<9l)] ■ 



(25) 



Since the operators I and J^, k = 2, . . . , /, mutually 
commute the eigenstates of Hqn F can be chosen such 
that they are simultaneously the eigenstates of all the 
elementary operators I and Jk , whose spectral properties 
are well known. Thus, 



HQ$ ¥ \I,n-z,...,nf) = E\I,n 2 , . . . ,n f ) (26) 



with 



\I,n 2 ,...,n f ) = \ijJi)®\ipn 2 )®---®\?Pn f ), (27) 
where 

J|Vi) = i"hM, (28a) 

J k \i> nk ) = h(n k + 1/2)| J , n k G N , (28b) 



and 



E = KqnA 1 ' Km + V2), ■■■,h(n f + 1/2)) . (29) 

Using the basis given by the eigenstates \I, n 2 , . . . , n/) 
one can now straightforwardly trace out the bath modes 
in Eq. (f2~4"|) . Indeed, let us define the operator 

[N/2] 

y(N) .fat 
n m,...,n f ~ n ' a '3 

n=0 \a\+j=n 

Q2 



712 + \ 



n/ + ' K\«\-°*+* (30) 



parametrized by the (/ — 1) nonnegative quantum num- 
bers n 2 , ■ • ■ , nf . Then Eqs (|24|) and (|25|) can be written 

as 

C QM {E, t)= Tr { S ( E - U)Fe iiit/h J : e~ i ' fLt/H } 

(31) 



«2, ■■■,tlf 



and 



(32) 



respectively, where, to avoid a cumbersome notation, we 
have dropped the superscript N and subscripts n 2 , . ■ ■ , rif 
for the operators H and T. 

Equations (|3"Tj) and (|3"2")l show that the problem of cal- 
culating the quantum FFCF for a system with / > 1 
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degrees of freedom effectively reduces to the correspond- 
ing problem for a one-dimensional system described by a 
Hamiltonian of the form Eq. (j30f which is a polynomial of 
the operator / associated with the reactive mode only. In 
the following section we present an explicit calculation of 
the FFCF for the simplest anharmonic one-dimensional 
Hamiltonian of this form. 



where for the discussion below, we explicitly added the 
anharmonicity parameter a to the argument of the FFCF 
(note that as opposed to a the parameter A can in prin- 
ciple be removed by a suitable scaling of the energy). If 
we denote the two solutions of E + aE 2 = E by 

1 



E ° = Y a (- 1 + °"v / l + 4a£) , tr = ±l, (39) 



III. THE QUANTUM FLUX-FLUX 
CORRELATION FUNCTION FOR ONE 
DIMENSIONAL ANHARMONIC BARRIERS 

We now present an analytical calculation of the FFCF 
for the two-parameter family of Hamilton operators de- 
fined by 

H(a,X) = ~h + ah 2 , h = - (f - X 2 q 2 ) . (33) 

Here A parametrizes the width of the barrier in the 
harmonic approximation and a characterizes the anhar- 
monicity of the barrier. The Hamiltonian operator H can 
be viewed to be in quantum normal form. In fact the op- 
erator h differs from the operator I defined in Sec. Ill Bl 
only by a factor of A which follows from a linear trans- 
formation of p and q which does not alter the normal 
form procedure described in the previous section. The 
Hamilton operator in (|3"3")l can therefore be considered to 
describe the simplest possible anharmonic barrier. 

The starting point of our calculation is the system of 
eigenstates of h, (and therefore of H) 



a = ±l 



(34) 



The corresponding wavefunctions are [20L l2ll 




. E 

'xh 



(35) 



where D v denotes the parabolic cylinder function of order 
v [22j ■ The eigenstates are mutually orthogonal, 



{i>^?,) = 8 rT<a ,8{E-E'), 
and form a complete basis, i.e., 



<T = ±1 



-foo 



dE\i^)(^\ = l. 



(36) 



(37) 



where 1 denotes the identity operator. 

Using the ip^ basis to expand the trace in Eq. (|3"Tj) we 
can write the quantum FFCF as 



C QM {E, t,a) = J J dE'dE" 8 (E" + aE" 2 - E] 



x exp 



E' + aE 12 - E) 



(38) 



then ([38]) becomes 
CQu(E,t, a) 



y/TTTa~E 7_ 



dE' e i(E'+aE' 2 -E)t/h 

(40) 



if f + AaE > 0, and CQu{E,t, a) = else. The latter 
condition on the energy E and the parameter a simply 
assures that we actually have a barrier scattering problem 
if the inequality is satisfied. 

The matrix elements of the flux operator F are calcu- 
lated as follows. According to Eq. (|32p we have 



F=-[h + ah 2 ,Q(q)] 



(41) 



so that 



(tf\F\1>$) =- [(E + aE 2 ) (E' + aE' 2 )] 
dq (^\q)(q\^'). 



(42) 



Then, using Eq. (|35|) into Eq. (j42|) and performing the 
integration over q we obtain 



x < ere 



87T 2 ft 

-**rf -h 



l + a(E + E') -K E + E> j E-E 

. E \ „ ( 1 . E' 



2Aft/ r U l 2XH 





. E > 




.E'\ 


a- 


Vl 2Xhj 


<\ 


l 2Xh 



^ r • (43) 



This finally leads to the following expression for the dou- 
ble sum (over a and a') entering Eqs. (|3"5|) and (jUI)) : 



[l + a(E + E')} 2 /nE + E' 

exp ( 2 ~Jn~ 



3 E 
l 4 +l 2Xh 



3 E' 
{ 4 +l 2Xh 



,1 E' \ 



(44) 



We then substitute Eq. ((33]) into Eq. ([3DJ and use the 
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formula (see Appendix [X] for its derivation) 

p+oo 

dxe lAx2+Bx {l + Cx) 2 \T(D + ix)\ 2 



n r ( 2£) ) 

= 2?r „„„ exp 



A d 2 



2 2D "4Qx 2 , 
(cosh x + iCD sinh x) 2 + \C 2 D 



(cosh x 



215+2 



(45) 



with A = AaX 2 ht, B = ir + 2iXt, C = 2Ata/(l + aE), 
and D = l/4or3/4to arrive at the central result of our 
paper: 



Cqm(E, t, a) = 29 X , e 2 J /2h HE/Xh, At, aXh) 



(46a) 



where 
A(e,T,a) 



VTTTe 



E 

<T=±1 



exp 




[ 3 - + l 6 - 

4 2 



-H 



I -;e 



(46b) 



e CT (e, a) = — (-1 + crVl + 4ea) , (46c) 
za 



and 

CT (^;£,r,a) 

2 ( d 2 \ (e a sinh r — 2iv cosh t) 2 — 2v 

= a J (sinhr) 2 ^ 2 • 

(46d) 

The expression for the FFCF given by Eq. (|4"6")l is ezact 
and holds for all energies E and parameters a satisfying 
1 + 4aE > 0. (As shown above, CQyi{E,t,a) = if 
1 + 4aE < 0.) In the following we consider the limits of 
an harmonic saddle, a — > 0, and short times At <C 1 for 
which cases VL a in Eq. (|46d|) assumes a simpler form. 



A. The case of a harmonic barrier (a — 0) 

In the limit a — > Eq. (|46cj) yields ei = s + 0(a) and 
e_i = — l/a — e + O(a). It is then straightforward to 
show that for a = the a = — 1 contribution to the sum 
in the right-hand side of Eq. (|46bj) vanishes and thus 



A(e,T,0) = e ; 



Ir (3 + is)| 2 Ir (i + ^ |2n 



(sinhr) 1 / 2 



(sinh r) 3 / 2 



(47) 
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FIG. 1: (Color online) The time dependence of the real (top) 
and imaginary (bottom) parts of the function A in (|47[) which 
enters the quantum FFCF for an harmonic barrier in (|46a[) . 
The graphs are for different values of the scaled energy e. 



which is valid for all energies, e € IR. 

Equations (|46a|) and (j47]) yield an exact expression for 
the microcanonical quantum FFCF of the parabolic bar- 
rier system with Hamiltonian (p 2 — X 2 q 2 )/2. The corre- 
sponding canonical version of the FFCF, defined as 



C QM (l3,t,a) = {e-Pn Fe iAt/h Fe- i6t/n ^ 



(48) 



can be readily calculated for the case of a = by perform- 
ing the bilateral Laplace transformation of (|46a|) which 
gives 

/+oo 
dEe- 0E C QM (E,t,O) 
-DO 

A 2 cosh(At c )sinh(Aft^/2) 



An [sinh 2 (At c ) + cosh 2 (A/l/3/2)] 



3/2 



(49) 



with t c = t- Here Eq. (g5j with A = C = was 

used to calculate the integral over energy. 

We note that Eq. (|4"9"|) was originally obtained by Miller 
et al. [l[ by representing the FFCF in terms of the time 
evolution operator for the harmonic oscillator. However, 
to our knowledge, the explicit expression for the micro- 
canonical FFCF in Eqs. (|46a|) and (|47[) has not been re- 
ported in the literature before. 

In Fig. m the time dependence of the function A in (|4"7) 
is shown for different values of the scaled energy e. We 
note that A and hence the microcanonical FFCF always 
diverge in the limit t — > and depending on the energy 
shows more or less pronounced oscillations. 
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FIG. 2: (Color online) Short time regime of the real (top) and 
imaginary (bottom) parts of the function A which enters the 
quantum FFCF in (|46a) (see Sec. HlTBj) . The graphs are for 
different values of the scaled anharmonicity parameter a and 
all have energy zero. 



B. The short-time regime (\t <C 1) 

For short times, r<l, one can approximate the hy- 
perbolic functions on the right-hand side of Eq. (|46d| by 
their leading order Taylor expansions to obtain 

VL a {v;e,T,a) 

~ a 2 (e 2 a V 2u - 4we a V 2l/+1 + 2v{2v + \)V 2v+2 ) (50) 



with 



d 2 \ i 



V„(t, a) = exp (-iar-^ j — . (51) 
As shown in Appendix [B] 2? m in Eq. (|5T|) can be written 



,2a 



1/2 



(52) 

Equations (|46aM6cj) together with Eqs. ([50]) and ([52]) 
provide an explicit expression for the FFCF at short 
times. 

Figure [2] compares the time decay of the dimension- 
less correlation function A in the harmonic case (a = 0, 
blue lines), given by Eq. (|47|l. and that in the anhar- 
monic case (a = 10 -3 , red lines), given by Eqs. (UObl) . 
([50]) , and ([52]) . One sees that even for very small (but 
non- vanishing) values of the dimensionless anharmonic- 
ity parameter a the time dependence of the FFCF signif- 
icantly differs from that of the corresponding harmonic 
problem. In fact, it is straightforward to show that as 
r — >• one has Re A —> +oo for a = 0, while RcA^ — oo 
for a > 0. The transition from the anharmonic case to 
the harmonic one takes place in a discontinuous manner: 



as a tends to zero the maximum of Re A (i.e., the peak 
of the red curve in the upper half of Fig. [3]) becomes 
higher and sharper and approaches r = recovering the 
harmonic result (the monotonic blue curve in Fig. [2]). 



IV. CONCLUSIONS 

In this paper we have presented a method for comput- 
ing classical and quantum flux-flux correlation functions 
for reactive systems with a potential barrier characterized 
by saddle type equilibria in phase space. The method is 
based on the normal form transformation of the system's 
Hamiltonian in a vicinity of the saddle point. 

In the classical case, the time dependence of the cor- 
relation function (with respect to a recrossing free divid- 
ing surface) is given by the (5-function. While this form 
of the time dependence has been known in the theoreti- 
cal chemistry community for some time, the contribution 
of the classical normal form theory is that it provides a 
dividing surface having the no-recrossing property that 
allows the computation of the pre-factor (essentially the 
flux through the dividing surface). No computation of 
trajectories is required to evaluate the flux-flux correla- 
tion function. The time integration to compute the rate 
becomes trivial. 

In the quantum case, we showed that the problem of 
calculating the correlation function in a system with more 
than one degrce-of-freedom reduces to an effective one 
degree-of-freedom problem. The Hamiltonian of this ef- 
fective one degree-of-freedom system is obtained through 
the quantum normal form procedure. Finally, and most 
importantly, we derive (for the first time in the literature) 
an analytical expression for the flux-flux correlation func- 
tion for the simplest anharmonic one-dimensional Hamil- 
tonian in quantum normal form. 
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Appendix A: Derivation of Eq. (|45[) 

We begin our derivation of the formula, Eq. (|45p . by 
writing 



ir(^>r=jp 



dq ■ 



sh(g/2)] 



2D 



(Al) 



The integral representation given by Eq. (1A1|) is readily 
obtained, e.g., from formula 5.13.2 in Ref. [2J]. Then, 



8 



denoting the left-hand side of Eq. (|45j) by 2" we get 
r(2D) [ + °° dq 



2 2D 



[cosh(g/2)] 



2D 



+00 



x / dx{l + Cx) 2 e iAx -<i+* B ) x . (A2) 



The second integral in the right-hand side of Eq. (|A2[) 
can be written as 

[ + °° dxx 2n {1+Cx) 2 e - l{q+lB )^ (A3) 



n=0 
Now, using 

"+00 



dxx n e~ t{q+tB)x = 2m n 6^ (q + iB) , (A4) 



with S( n ' denoting the nth derivative of the delta func- 
tion, and then, carrying out the g-intcgration in Eq. (|A3[) 
we obtain 

_ _ r(2,D) ~ d 2n 
1 = 2tt^^- ^ 



2 2D ^ n\ dq 2n 

n=0 ^ 



coshf -MCT>sinhf) + \C 2 D 



, 2D+2 



(coshf) 

Finally, formally summing the series 
(-iA) n d 2n 



(A5) 



q=-iB 



£ ^if^^ f{q) = cxp I l 



n=0 



with / denoting an arbitrary function, and making the 
change x = q/2 we arrive at Eq. (|45|) . 



Appendix B: Derivation of Eq. (|52[| 



Taking into account the identity 



1 



we rewrite Eq. (|5ip as 



7? — 

" r( M ) 



Then, using formula 3.462.1 in Ref. [24J, 



r( u ) / b 2 , 



B 

7Tb 



(Bl) 



— / dss^e- 101 ™ 2 -™ . (B2) 
^) Jo 



(B3) 



for A = e"/ 2 aT and B = t we arrive at Eq. (|52l . 
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